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11 Abstract 

12 A neutron counter designed for assay of radioactive materials has been adapted for beam exper- 

13 iments at TUNL. The cylindrical geometry and 60% maximum efficiency make it well suited for 

14 (7, n) cross-section measurements near the neutron emission threshold. A high precision charac- 

15 terization of the counter has been made using neutrons from several sources. Using a combination 

16 of measurements and simulations, the absolute detection efficiency of the neutron counter was de- 

17 termined to an accuracy of ± 3% in the neutron energy range between 0.1 and 1 MeV. It is shown 
is that this efficiency characterization is generally valid for a wide range of targets. 
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20 I. INTRODUCTION: 



21 The model IV inventory sample counter (INVS) developed at Los Alamos National Lab- 

22 oratory [1J was designed for fast, non-destructive assay of radioactive materials. Specialized 

23 inserts for the axial bore of this neutron counter have been made to adapt it for use as 

24 the primary neutron detector for in-beam (7,n) total cross section measurements. Develop- 

25 ment and testing of this counter took place at TUNL using hadron beams in the tandem 

26 laboratory and the 7-ray beam at the High Intensity Gamma-ray Source (HI7S) [2]. 

27 Such measurements require detailed and accurate information about the energy-dependent 

28 absolute neutron detection efficiency of the counter. Efficiency here is generally defined as 

detected / \ 

€ = N ' (1) 

1 y emitted 

30 The efficiency measurements were made using four different sources, each with a precisely 

31 known neutron emission rate. First, a 252 Cf source, calibrated by the National Institute 

32 of Standards and Technology (NIST), generated a flux of neutrons known to ± 4.4% [3]. 

33 Second, a coincidence experiment using the 2 H(<i, n) 3 He reaction provided a mono-energetic 

34 source of 2.26 MeV neutrons with flux known to ± 10%, and gave insight into the thermal- 

35 ization time of neutrons in the INVS. Third, an investigation of the 7 Li(p, n) 7 Be reaction 

36 produced <1 MeV neutron sources with fluxes known to ± 6.6%. Finally, the 2 H(7,n) 1 H 

37 reaction was used to produce tunable sources of monoenergetic neutrons (0.1 < E n < 1.0 

38 MeV) with fluxes known to ± 3% accuracy. A comparison of all experimental data with 

39 simulations demonstrates varying levels of agreement. 

40 The following sections describe the detector geometry (Sect. [IT]) followed by details for 



each experimental setup (Sect. Ill) including discussions of backgrounds, calculations, mea- 



42 surement uncertainties and results. Section IIVI contains a discussion of Monte Carlo simu- 

43 lations used for comparison with experimental results. Finally, Sect. [V] offers a summary 

44 of the results and a discussion of future applications of the INVS counter for cross-section 

45 measurements. 



46 II. DETECTOR DESCRIPTION: 



47 

48 



The active detection elements in the INVS counter are 18 tubular proportional counters, 
each containing 6 atm. of 3 He. The tubes are arranged in two concentric rings at radii 7.24 
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49 cm and 10.60 cm each containing nine equally spaced detectors (see Fig. [T]). The detectors 

50 are embedded in a cylindrical polyethylene body 46.2 cm long and 30.5 cm in diameter 

51 which serves as a neutron moderator. The active length of the 3 He gas within the tubes 

52 is 39.4 cm pp. The detector body has an 8.9 cm diameter axial cavity designed to contain 

53 the neutron source. Throughout this manuscript the term longitudinal center refers to the 

54 center of the detector with respect to the length of the detector body, and is distinguished 

55 from the term axial center which refers to the axis of the detector. 

56 In experiments with beam, the irradiated target is the source of neutrons. During all but 

57 one of the experiments described here, the neutron-emitting source was located inside the 

58 detector cavity, usually near the longitudinal center. The cavity was partially filled with 

59 additional neutron moderator (often graphite and/or polyethylene) to increase the detection 
eo efficiency. The additional moderator was arranged so that the beam could pass through the 
ei detector without intercepting moderator material. 

62 Thermalization of the neutrons within the detector body increases the probability for 

63 initiating the 3 He(n,p) 3 H reaction within the embedded tubes. An energy of 763.7 keV, 

64 shared between the outgoing proton and triton, is released from each reaction. Most of the 

65 kinetic energy is lost to ionization of 3 He, which is detected as an electrical pulse on the 
ee central electrode of each tube which is biased to +1780 V. A fixed threshold effectively dis- 

67 criminates against low-pulse-height signals generated by 7-rays and electronic noise. Signals 

68 above the threshold generate ~50 ns wide TTL pulses using on-board electronics. Onboard 

69 signal-processing electronics within the detector produce three TTL logic output signals; the 

70 inner ring (/); the outer ring (O); and the logical OR of the I and O pulses (T). Whenever 

71 one or more tubes in the inner (outer) ring detect a neutron, a pulse is generated on the / 

72 (O) output. For neutrons with energies less than about 2 MeV, the I/O ratio can provide a 

73 reasonable determination of the mean of the energy distribution of the detected neutrons. 

74 III. EXPERIMENTS: 

75 In this section, the experimental setup and techniques used to measure the efficiency of 

76 the counter are described, and the results of the each measurement are presented. Discussion 

77 of the results will be presented in the next section. Four different sources of neutrons were 

78 used to cover an energy range from about 0.1 to 10 MeV. 
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FIG. 1. (Color Online) Front and side cut-away cross-sectional views of the Model IV INVS counter 
described in the text. The arrangement of inter-cavity moderator corresponds to the experimental 
geometries either for the 252 Cf source measurement, or for the 2 H(7, n) 1 !! experiment. 
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FIG. 2. (Color Online) Efficiency vs. Z-axis position for an open detector geometry. 
A. 252 Cf 

so Californium-252 is a standard calibration source for neutron detectors. The effective half 
si life of 252 Cf is 2.645 years and is due to alpha particle emission and spontaneous fission which 
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82 produces a neutron yield of 2.314 x 10 6 neutrons/s//zg [Ij. The energy spectrum of neutrons 

83 from 252 Cf is well known [5]. A calibrated 252 Cf source provides a single measurement of 

84 efficiency representing the response of the detector to a broad spectrum of neutron energies. 

85 The 252 Cf source used consisted of 3.15 ng of active powdered material encapsulated 

86 within a small aluminum pellet. The source was suspended on the axis of the detector 

87 approximately 3.4 cm from longitudinal center. Graphite moderator filled most of the volume 
as of the cavity. A table of neutron activity versus date was provided by NIST [3]. The 

89 experimentally determined efficiency for this configuration (shown in Fig. [TJ was 40.5 ± 1.8 

90 %. The experimental I/O ratio was 1.516 ± 0.004. 

91 The dependence of the detection efficiency on the position of the source within the central 

92 cavity was determined by making measurements with the source placed at different positions 

93 within the central bore. Measurements on the central axis were made along the entire length 

94 of the detector. The detection efficiency has a maximum value at the longitudinal center and 

95 drops off smoothly as the source is moved in either direction away from the center along the 

96 detector axis. The shape of the position dependency is a purely geometric acceptance effect 

97 and can be approximated analytically for point sources with isotropic neutron emission. The 

98 measured detection efficiency as a function of the source position along the central axis of 

99 the counter is shown in Fig. [2] and in comparison to simulated and calculated efficiencies. 

100 The relative efficiency is directly proportional to the angular acceptance of the counter 

101 as a function of z, which is given by the equation below for an isotropic point source of 

102 neutrons. 




L/2 



v/(L/2 - zf + r 2 

L/2 + z 



+ 1 



(2) 



^{L/2 + zf + r 

we Here, r is the radius of the opening at the end of the detector, L is the active length of 
107 the 3 He gas, and z = at L/2. This function is maximum when z = 0. For an open cavity 
los geometry, the change in the efficiency over the length of an 8-cm long sample, centered 
109 on the axis at the longitudinal center is approximately 1%. For a geometry like the one 
no shown in Fig. [T] the changes in efficiency are negilgibly small over a length of nearly 20 cm 



in centered on the longitudinal center. For sources located off its central axis, the detection 

112 efficiency changes by less than 0.5% . The radial dependence of the efficiency is also mostly 

113 a geometric acceptance effect. 



ii4 B. 2 H(d,n) 3 He 

H5 The 2 H(d, n) 3 He reaction was used to measure the efficiency for monoenergetic 2.26 MeV 

lie neutrons. The associated particle technique was used with the recoil 3 He nucleus detected 

117 in a siliFcon surface barrier detector inside an evacuated chamber. A schematic diagram of 

us the experiment setup is shown in Fig. [3j The neutron counter was positioned so that its 

n9 central axis coincided with the symmetry axis of the cone of neutrons associated with the 

120 3 He particles detected in the silicon detector on the opposite side of the incident beam axis. 

121 The distance from the longitudinal center of the counter to the deuterium target was set so 

122 that the diameter of this neutron cone was smaller than the diameter of the central cavity 

123 through the detector. The energy of the incident deuteron beam and the detection angle of 

124 the silicon detector were set to produce 2.26-MeV neutrons emitted along the central axis 

125 of the counter. With this method the efficiency is computed as 

(3) 



Ns He 

127 where N3# e is the total number of detected 3 He-particles and N n is the total number of 

128 neutrons detected in coincidence with the detected 3 He particles. This equation takes the 

129 detection efficiency of the silicon detector to be unity. The deuterium targets used in these 

130 measurements were ~ 100/ig/cm 2 thick deuterated polyethylene (C2D4) evaporated onto a 

131 10 /ig/cm 2 thick carbon foil. The deuteron beam energy incident on the foil was 2.0 MeV, 

132 and the average beam current on the C2D4 foil was ~20 nA. The cross-sectional profile of 

133 the deuteron beam at the foil was circular with a diameter of approximately 0.5 cm. Each of 

134 the two silicon detectors (one in-plane and one out-of-plane) had a solid angle acceptance dfl 

135 = ?t/60 sr, and each was located at a scattering angle of 9^ = 26.50°. Neutrons associated 

136 with detection of 3 He in the in-plane Si detector exited the target at 6^ = 117.1° along 

137 the central axis of the neutron counter. The rear half of the central cavity was plugged 

138 with polyethylene to scatter neutrons traveling through the central cavity into the body 

139 of the counter. The out-of-plane Si detector was used to measure the rate of accidental 



Si Detectors 



FIG. 3. (Color Online) Schematic diagram of the experimental setup for the efficiency measure- 
ments made using the 2 H(d, n) 3 He reaction at the tandem accelerator facility. 
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FIG. 4. (Color Online) The histogram is a background subtracted TAC spectrum from the 
2 H(d, n) 3 He experiment. A simulated TAC spectrum (blue points) is shown in comparison to 
the experimental histogram. The solid curve is a fit to the simulated spectrum and predicts that 
most neutron detection occurs outside the 22.5 /xs experimental window. 

wo coincidences. 

141 Efficiency and thermalization time were deduced simultaneously using a time-to-amplitude 

142 converter (TAC) which recorded the time between a charged particle detection in the silicon 

143 detector and a neutron detection in the IN VS. A threshold setting effectively discriminated 

144 against deuteron elastic scattering events. The effective TAC range was 22.5 /is. The 

145 TAC was calibrated using a pulser which started and stopped the TAC with known delay, 
we Because the INVS is a thermalization counter, detection efficiency is time-dependent on a 

147 microsecond time-scale. A peak in the TAC spectrum at ~3 /xs suggests a source of delay 

148 exists caused by charge collection and signal processing in the INVS counter. This delay 

149 reduces the effective TAC range to 22.5 /zs from 25.5 /zs, which is where the experimental 
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150 TAC spectrum ends. The present result for neutron detection efficiency over a 22.5 /is range 

151 is 11.0 ± 1.1%. 



152 C. 7 Li(p,n) 7 Be 

153 The 7 Li(p, n) 7 Be reaction was used to measure the energy-dependent detection efficiency 

154 over an energy range that overlaps with that covered by the 2 H(7, n) 1 !! source reaction below 

155 about 0.7 MeV and to provide data for a neutron source with the intensity distribution 

156 peaked at forward angles relative to the central detector axis P, |7j. The cross section for 

157 7 Li(jo, n) 7 Be reaction is large and has been accurately measured [8] making it a good neutron 

158 source for calibrating the efficiency of detectors at low energies [[JJ. 

159 1. Experimental Setup 

wo The experimental arrangement is shown in Fig. [5] The proton beam was tuned through 

lei a double collimator set onto the LiF neutron production target. The cross-sectional profile 

162 of the beam on target was circular with a diameter of 5 mm, and the average beam current 

163 on target was 100 nA. The energies of the proton beams incident on the LiF target were 

164 between 1.88 and 2.46 MeV. The neutron production target was comprised of 39.8 /ig/cm 2 

165 of LiF evaporated onto a 8.3 /zg/cm 2 thick carbon backing. Targets were located on the axis 
we of the INVS counter inside an evacuated beam pipe at 14.2 cm from the longitudunal center. 
167 The transmitted proton beam was collected in a voltage-suppressed Faraday cup at the end 
lea of the beam pipe. A polyethylene plug was placed just beyond the end of the beam pipe to 
169 increase detection efficiency. Backgrounds were measured by putting beam through both an 
no empty target ring identical to the one that supported the LiF target, and a target ring that 

171 supported only a carbon backing. In total, beam-induced and environmental backgrounds 

172 amounted to < 0.1% of real counts. 

173 2. Results 

174 The detection efficiency as a function of proton energy was calculated using 
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Collimator Faraday Cup 

FIG. 5. (Color Online) Schematic diagram of the experimental setup for the efficiency measure- 
ments made using the 7 Li(p, n) 7 Be reaction at the tandem accelerator facility. 
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FIG. 6. (Color Online) For the 7 Li(p, n) 7 Be reaction, a comparison of e(-Ep) as determined by 



experiment (points) and by simulation (band). Simulations, discussed in Sect. IV, reproduce the 
shape of the efficiency well. 



176 where N n is the total number of neutrons detected, N p is the number of protons collected 

177 in the Faraday cup, N t is the number of target nuclei per unit area, and <j(E p ) is the total 

178 cross-section of the 7 Li(p, n) 7 Be reaction at proton energy E p . 

179 The data (see Fig. [6]) display a relative minimum in efficiency near E p = 2.13 MeV followed 
wo by a relative maximum near E p = 2.32 MeV. These shifts in efficiency coincide with rapid 

181 changes in the angular distribution of neutrons. Though statistical uncertainties were very 

182 small, systematic uncertainties for target thickness and cross-section contributed 3.5% and 

183 5%, respectively, resulting in an overall systematic uncertainty of 6.6%. 
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FIG. 7. (Color Online) Schematic of the experimental setup for measurements made using the 
2 H(7,n) 1 H reaction at the HI7S facility. After collimation, the 7-ray beam passes through scin- 
tillation paddles (not shown) and into the target room. The 1.2 cm diameter 7-ray beam passes 
through the following elements: (a) "clean-up" collimator wall; (b) the D2O target located near 
the longitudinal center of the neutron counter; (c) machined lead attenuators located between lead 
collimator walls; (d) a Nal detector; (e) an HPGe detector. 

184 D. 2 H(7,n) 1 H 

185 The 2 H(7,n) 1 H measurement was unique among the experiments described here in that 
we it produced nearly monoenergetic neutrons with very small flux uncertainties. Several ef- 

187 ficiency measurements were made that highlighted the energy-dependent response of the 

188 detector. A schematic of the experimental setup is shown in Fig [7} 

189 1. Experimental Setup 

190 The target consisted of approximately 3.2 g of 99.8% enriched D 2 sealed inside a thin- 

191 walled polyethylene tube that was 7.62 cm long. The target was located at longitudinal 

192 center, but 2.9 cm above the axial center, and occupied one of four 2.54 cm diameter az- 

193 imuthal holes in a rotatable graphite cylinder. The other holes were available for other 

194 targets. A graphite target and an empty hole were alternately rotated into the beam during 

195 experiments to determine beam-induced backgrounds. The samples were remotely rotated 

196 into position using a four-position Geneva mechanism which assured reproducible alignment 

197 on the 7-ray beam axis. Most of the data were collected for 7-ray beams with energies 

198 between 2.48 and 4.10 MeV. 
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199 Gamma-ray production by inverse- Compton scattering at HI7S is well documented [2, 

200 Q2I1 E]. The 7-ray beams used in this experiment were collimated to 1.2 cm in diameter. 

201 The size of the beam, and its alignment with the target was confirmed using a 7-ray beam 

202 imaging system [12] . 

203 The 7-ray beam energy distribution was determined by a high purity germanium (HPGe) 

204 detector located on the 7-ray beam axis. A radioactive 60 Co source, and naturally present 

205 40 K and 208 T1 provided energy calibration. The FWHM of the 7-ray beam as determined 

206 by the HPGe was typically between 1 - 3%. A typical 7-ray spectrum is shown in Fig. |8j 

207 Relative incident 7-ray flux was continuously monitored by three scintillating paddles 

208 located upstream from the experimental setup. The absolute 7-ray fluxes were determined 

209 by a cylindrical 25.4 cm x 35.6 cm Nal detector located behind the active target, on the 

210 7-ray beam axis. For the 7-ray beam energies for which it was used, the Nal detector had 

211 a total integrated efficiency of nearly 100%; In other words, nearly all 7-rays will interact 

212 within the detector volume and deposit energy. A threshold setting which ignored signals 

213 generated from 7-rays with energy less than ~0.6 MeV reduced the efficiency to ~97% for all 

214 experimental 7-ray energies as determined by simulations. This threshold setting optimized 

215 the ratio of total integrated signals to room background. 

216 After passing through the target, the 7-ray flux was attenuated by machined lead atten- 

217 uators to eliminate signal pile-up and dead-time effects. This allowed 7-ray flux on target 

218 of >10 6 7/s, and flux at the face of the detector of <10 4 7/s. Determination of the absolute 

219 beam flux from the signals measured in the Nal detector requires precise quantitative infor- 

220 mation about the effective attenuation of the beam by the lead attenuators at each 7-ray 

221 energy. The effective 7-ray attenuation, was determined for each attenuator, at several 7-ray 

222 energies, using the scintillating paddle system for flux normalization. 

223 2. Background 

224 Two sources of neutron backgrounds existed in this experiment. One source is classified as 

225 environmental backgrounds. These can be caused by either cosmic-ray neutron production 

226 or natural radioactive sources in the vicinity of the INVS counter. These sources generated 

227 0.2 count/s per 3 He tube for a total of 3.6 counts/s. The second background source arose 

228 from 7-ray-beam-induced events. Gamma-rays that scatter from the target can deposit 
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229 enough energy to register a signal above the threshold setting. This type of background 

230 was measured by bombarding a graphite target with the 7-ray beam. Gamma-ray beam- 

231 induced backgrounds amounted to approximately 4.8 counts/10 6 7-rays on target, which 

232 was typically < 1% of real counts from the heavy water target, and were taken into account 

233 in the data analysis. 



234 3. Results 



235 For the experimental setup used in these measurements the detection efficiency for neu- 

236 trons from 7-rays on a heavy water sample can be explicitly calculated from 

{E) = N nX { El )e{E) 

K <' fN^N t a{Ej 1 } 

238 where N n is the number of neutrons detected, x(E 7 ) is the measured 7-ray attenuation by 

239 the lead attenuator at 7-ray energy E 7 , e 7 ( E 7 ) is the efficiency of the Nal detector for 

240 7-rays with energy E 7 , f is the thick target correction factor (described below), iV 7 is the 

241 number of 7-rays detected, N t is the number of target nuclei per unit area and cr(E^) is the 

242 total cross-section of the 2 H(7,n) 1 H reaction at E 7 [13]. The author of Ref. [13] calculated 

243 the 2 H(7,n) 1 H total cross-sections with several widely-used N-N potential models, all of 

244 which were indistinguishable to within 1%, irrespective of the model used. In addition, this 

245 cross-section agrees with the world data which have uncertainties between 3 and 6% [HHT7] . 

248 

248 The f factor in Eqn. [5] accounts for the flux loss due to interactions with atomic electrons 

249 as the 7-rays propagate through the heavy water target. This factor is calculated as, 

1 \ f(l - e-"»*) 



251 using NIST attenuation coefficients for heavy water /x w [18] and the mass thickness t of the 

252 target. 

253 For photodisintegration of the deuteron E 1 and E n are related, in units of MeV, by 

E., - 2.225 

2.001 V ; 

255 where -2.225 is the Q-value for the reaction, and the factor of 2.001 comes from energy 

256 sharing between the outgoing proton and neutron. Cross-sections for energies 2.48 MeV < 

257 E 7 < 4.10 MeV were used to ensure ± 1% accuracy. 
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FIG. 8. (Color Online) A HPGe spectrum for 2.470 MeV 7-ray beam with dE/E = 1%. Spectra 
of calibration 7-rays from 60 Co, 40 K, and 208 T1 are overlaid. 

258 The total combined statistical uncertainty in the efficiency measurements made using 

259 neutrons from deuteron photodisintegration was < 2%. The main source of statistical error 

260 was in the attentuation measurements of the lead absorbers. The statistical accuracy of 

261 these measurements was determined by the counts in the scintillator paddle system used 

262 for relative flux normalization between the absorbers. A minimum of ± 1% statistical 

263 uncertainty was achieved for most energies. The statistical uncertainties in the 7-ray beam 

264 flux measurements and the neutron counting for the heavy water target were typically < 

265 0.5%. 

266 The total systematic uncertainty was < 3% and was mainly due to three sources. The 

267 uncertainty in the deuteron photodisintegration cross section was kept to less than 1% by 

268 limiting the energy range of the measurements to 2.48 MeV < E 7 < 4.14 MeV. Uncertainties 

269 in the target thickness and cross section contribute 0.5% and 1% respectively. 




270 IV. DISCUSSION AND RESULTS: 

271 A. Simulations: 

272 The Monte-Carlo code mcnpx was used to simulate all particle interactions within the 

273 INVS counter for each experiment. For all simulations, material densities for the INVS 

274 counter were fixed, and standard thermal neutron capture cross-section libraries (19] were 
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275 used. Variable parameters in each simulation were: 

276 1. the arrangement of materials inside the cavity; 

277 2. the location of the neutron emitting source; 

278 3. the energy and spatial distribution of neutrons. 

279 Absolute detection efficiencies and I/O ratios were extracted from simulations for comparison 

280 with each experiment where applicable. 



252 



Cf 



282 Efficiency measurements of 252 Cf were made during the course of the 2 B.( / ~f,n) 1 B. exper- 

283 iment. Consequently, the arrangement of materials inside the INVS cavity was identical 

284 for both experiments. The energy distribution of neutrons produced by fissioning 252 Cf was 

285 modeled as a Watt fission spectrum which has the form 

p(E) = Cexp { - E/a) sinh{bE) 1/2 , (8) 



287 where a and b are parameters given for 252 Cf 

288 The MCNPX-simulated efficiency of 39.2% agrees with the experimentally determined 

289 efficiency of 40.5 ± 1.8 %. The simulated I/O ratio of 1.59 falls short of agreement with the 

290 experimentally determined 1.516 ± 0.004 due to a 6% larger efficiency for measurement in 



29i the outer ring (see Fig. 12). The source of this discrepancy is unclear. 



292 2. 2 H(d,n) 3 He 

293 The simulation of the measurements made with the 2 H(d, n) 3 He reaction were made as 

294 follows. A simulated beam of 2.26 MeV neutrons was emitted from inside an evacuated 

295 volume, through an aluminum beam-pipe wall, directed toward the axial center of a set of 

296 polyethylene plugs that filled most of the detector cavity (see Fig. [3]) . Because the experiment 

297 only recorded counts in the neutron detector during a 22.5 /xs wide time window after the 

298 associated 3 He particle was detected, it was necessary to track the neutron detection time 

299 in the simulations. Therefore, a time dependent model for detection of neutrons emitted 

300 from the 2 B.(d, n) 3 He reaction was created to compare with experiment. In this model, only 
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FIG. 9. (Color Online) (Top) A three parameter fit (solid line) describes the simulated efficiency 
of neutron detection as a function of time in the 2 H(c?, n) 3 He experiment very well. 
(Bottom) The same plot expanded to show t < 35 /xs. The vertical and horizontal lines identify 
the experimental window and expected efficiency 



301 neutrons detected before a user-defined time counted toward efficiency. A plot of efficiency 

302 vs. time was simulated for times between t = and t = 1000 /xs (see Fig. [9J. To produce a 

303 simulated TAC spectrum, a plot of the slope of e(t) vs. time was generated for comparison 

304 with data (see Fig. [i]). 

305 The simulated total efficiency for neutrons collected between and 22.5 /is is 11.9% in 

306 agreement with experiment. It is noteworthy, that simulations predict a relatively long time 

307 (nearly 500 /xs) before a maximum efficiency of 38.8% detection is realized for 2.26 MeV 

308 neutrons from the 2 H(<i, n) 3 He reaction. 
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3. 7 Li(p,n) 7 Be and 2 H(-/,n) l H 



310 Simulations for the 7 Li(p, n) 7 Be and 2 H(7, n) x H reactions were carried out in the following 

3n way. First, the location of the source was set to match experimental conditions. For a single 

312 simulation the source emitted monoenergetic neutrons only between angles 9 and 9 + dd with 

313 constant emission over <ft. After stepping through all of 9 space, the process was repeated 
3u for a new neutron energy. 

315 Ultimately, a three-dimensional plot was constructed with neutron energy on the x-axis, 

316 emission angle on the y-axis and detection efficiency on the z-axis (see Sect. [V]). After 

317 choosing an incident particle energy, and inputing expected angular distributions for the 

318 neutrons in the center-of-mass (CoM) frame, a second Monte Carlo process produced an 

319 average efficiency for the given source conditions. This process was repeated for several 

320 incident particle energies, and the result was a plot of simulated efficiency as a function of 

321 incident particle energy. 

323 For the 7 Li(p, n) 7 Be reaction, simulations reproduce very well the shapes of e(E p ) for both 



324 the inner and outer detector rings (see Fig. 11). Absolute detection efficiency for the outer 



325 ring of detectors is in good agreement with experiment. A 13% systematic offset in absolute 

326 efficiency is observed for the inner ring of detectors. Known systematic effects can account 

327 for a maximum difference of 6.6%. The difficulty with this discrepancy is that it appears 

328 to be of a systematic nature, while only affecting the inner ring of the INVS counter. A 

329 missing thermal neutron sink in the model may explain the difference. An underestimated 

330 amount of aluminum in the modeled beam pipe could have inadequately converted neutrons 

331 to 7-rays, leaving an excess of thermal neutrons in the region of the inner ring of the INVS 

332 counter. 

333 For the 2 H(7,n) 1 H experiment plots of e(E n ) for the inner ring, the outer ring, and the 

334 total show the data trends in good agreement with trends predicted by the simulation. The 

335 data for total detection efficiency are systematically lower than simulation by about 5.9%. 

336 Data for the inner and outer rings were systematically lower by 6.7% and 3.7% respectively. 

337 These data provide a benchmark calibration for this INVS counter with regards to its use 

338 in future (7, n) experiments. 

339 In the INVS counter a single detected neutron provides no information about the energy 

340 of the neutron. However, the average neutron energy from an ensemble of detected neutrons 
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FIG. 10. (Color Online) The efficiency (a) and I/O ratio (b) are shown for neutrons detected 
during the 2 H(7,n) 1 H experiment. Statistical error bars are smaller than the data points. 
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FIG. 11. (Color Online) Absolute efficiency for the inner ring (squares) and outer ring (triangles) 
of the detector during the 7 Li(p, n) 7 Be experiment. Statistical uncertainties are smaller than the 
data points. Simulations with uncertainties are shown as colored bands. 
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341 may be gleaned from the observed proportionality 

342 IfO(xE~K (9) 

343 which is easily inverted. The ability to distinguish the signature I/O ratio for E n from the 



344 I/O ratio for E n + AE n becomes more difficult as E n increases (see Fig. 10). 

345 The method described above for evaluating detector efficiency is valid irrespective of the 

346 target used. A deuteron target was chosen because of the precision with which the (7, n) 

347 cross section is known for this nucleus. Evaluating the efficiency this way established the 

348 energy-dependent response of the detector for future users under certain conditions. The 

349 efficiencies determined here are valid for any (7, n) reaction measurement that satisfies the 

350 following conditions: (a) the target location was fixed to match simulation; (b) the average 

351 energy of the emitted neutrons is known; and (c) the CoM angular distribution is known. If 

352 backgrounds are low, and the neutron energies are < 2.0 MeV, condition (b) may be relaxed, 

353 because the neutron energy information may be obtained from the I/O ratio. For neutron 

354 energies < 500 keV, condition (c) may also be relaxed because the response of the detector 

355 is nearly constant with respect to angle of emission. 



356 V. SUMMARY AND CONCLUSIONS: 

357 The goal of this work was to characterize precisely the response of the highly efficient 

358 INVS counter using multiple neutron sources with a focus on neutrons of energy < 1.0 MeV. 

359 The attention to low energy neutrons was motivated by a need to generate high quality 

360 (7, n) cross-section data. Experiments were carried out and then simulated in detail for 

361 comparison. In all simulations, special attention was given to assure accurate reproduction 

362 of experimental conditions. 

363 Figure [12] shows the ratio of experimentally determined efficiency to simulated efficiency 

364 vs. average neutron energy for each experiment. Ratios have been determined for the inner 

365 (I) and outer (O) rings separately, as well as for the total (T). 

366 Neutrons from the 2 H(7, n) 1 !! reaction were emitted from within ± 4 cm from the longi- 

367 tudinal center of the detector with a sin 2 (#) distribution in the CoM frame which is hardly 

368 changed when converted to the lab frame because of the relatively small momentum of 

369 the incident 7-ray. The arrangement of moderating materials was approximatly symmet- 
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FIG. 12. (Color Online) Plot of the ratio of data to simulation for all experiments. Results for 
total e n , inner ring e n , and outer ring e n are shown. Blue filled circles are 2 H(7,n) 1 H data; red 
squares are 7 Li(p, n) 7 Be data; black triangles are 252 Cf data; open circles are 2 H(d, n) 3 He data. 




FIG. 13. A three-dimensional plot of efficiency vs E n and Q[ a b using a likely arrangement of 
moderator for (7, n) experiments at HI7S. The simulated neutron-emitting target was at the axial 
and longitudinal center of the detector. The simulated geometry filled the detector cavity with 
graphite except for a 2.54 cm diameter hole for the target and 7-ray beam. 

370 ric. Under these conditions, simulations reproduce measurements to within a normalization 

371 of -6.7% (I), -3.7% (O), and -5.9% (T). These differences are likely the result of effec- 

372 tive threshold settings on the signals from the detector tubes. This feature of the effective 
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373 threshold settings was not included in the simulations. 

374 Neutrons from the 7 Li(p, n) 7 Be reaction were emitted far from the longitudinal center 

375 of the detector with a Legendre polynomial series distribution in the CoM frame which is 

376 significantly changed when converted to the lab frame because of the relatively large mo- 

377 mentum of the incident proton. The arrangement of moderating materials was asymmetric 

378 - biased to achieve higher detection efficiency for neutrons emitted near 9i a b close to zero. 

379 Under these conditions, absolute detection in O is reproduced very well by simulations; how- 

380 ever, absolute detection in I is systematically 13% smaller than predicted by simulations, 

381 resulting in a nearly 10% systematic difference in T. 

382 The difference in the level of agreement between simulations and experiment for I and O 

383 for the 7 Li(j>, n) 7 Be reaction is intriguing. The most likely explanation is that the amount of 

384 aluminum in the intervening beam pipe was underestimated in the simulations which caused 

385 an excess of thermal neutrons in the vicinity of /. 

386 The systematic differences between experiment and simulation for /, O and T for the 

387 2 H(7,n) 1 H measurement confirmed the need for a well known, tunable, monoenergetic neu- 

388 tron source. Reliance on simulations alone would have introduced systematic errors in future 

389 (7, n) measurements on the order of 6%. Using the cross-section of Ref. [13] as a 1% standard 

390 provided tunable monoenergetic neutron sources with fluxes known to ± 3%. The method 

391 for simulating absolute detection efficiencies for the 2 H(7, n) 1 H reaction was not significantly 

392 influenced by the choice of target material. Thus, using the same techniques, other (7, n) 

393 cross-sections may be measured with very high accuracy. 

394 For the purpose of planning future (7, n) experiments at HI7S which seek to use the INVS 

395 counter, simulations have been prepared for three likely experimental setups. An analysis 

396 program takes as input the geometry, the desired (7, n) reaction, the incident particle energy, 

397 and the CoM angular distribution of the outgoing neutrons and gives as output the expected 

398 detection efficiency and I/O ratio. The capabilities offered by the combination of the HI7S 

399 facility and this precisely characterized INVS counter make possible absolute photonuclear 

400 cross-section measurements with high precision. 
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